Sudden death of effective entanglement 
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Sudden death of entanglement is a well-known effect resulting from the finite volume of separable 
states. We study the case when the observer has a limited measurement capability and analyse 
the effective entanglement (i.e., entanglement minimized over the output data). We show that in 
the well defined system of two quantum dots monitored by single electron transistors, one may 
observe a sudden death of effective entanglement when real, physical entanglement is still alive. For 
certain measurement setups, this occurs even for initial states for which sudden death of physical 
entanglement is not possible at all. The principles of the analysis may be applied to other analogous 
scenarios, such as estimation of the parameters arising from quantum process tomography. 
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I. INTRODUCTION 

Sudden death of quantum entanglement [H-Q is one of 
the phenomena related to the fact that in finite dimen- 
sional systems the set of nonentangled states is of finite 
volume The phenomenon was explicitly identified in 
Ref. [l| J'see also its implicit presence in independent 
analysis Q). Full recognition of its iniportance and con- 
sequences was established with time [I i i H (for the 
review see Ref. @]) and demonstrated experimentally 0. 

In earlier analysis of entanglement evolution, it was 
assumed that the observer has the power to perform ar- 
bitrary measurements and can determine the state of the 
system completely. This is however not always true. In 
particular, there are natural systems, like quantum dots 
(QDs) which we will consider later, where limited mea- 
surement power is a natural and practical constraint (ie. 
single electron transistors (SETs) coupled to QD systems 
can be used to find a limited amount of information about 
the QD state [13, [13) In all such cases of limited measure- 
ment capability, it is natural to consider the worst case 
scenario proposed in Ref. [l2| : as real entanglement, one 
should consider the entanglement (i.e. chosen entangle- 
ment measure) minimized over the set of measurement 
data. This approach has found wide developments in 
terms of entanglement witnesses [l3l - [l6t with respect to 
experimental data [13, [13 • The minimized entanglement 
will subsequently be called effective entanglement. 

Here, we consider the system of a double QD (DQD) 
interacting with a phonon-bath which leads to an un- 
avoidable partial pure dephasing effect typically on pi- 
cosecond timescales [l9l - [21| . The time-evolution of phys- 
ical entanglement in this system can be found in Ref. 
[121 . We take into account the limitations to the knowl- 
edge of the system state imposed by a realistic measure- 
ment setup (consisting of different configurations of SETs 
interacting with the double-dot system). This measure- 



ment scheme does not allow for state tomography and 
provides, in fact, a very limited set of observables. In 
each time-step we minimize the value of entanglement 
with respect to the data which can be measured and 
find the evolution of the effective entanglement with re- 
spect to the SET-defined observables. We show that such 
an evolution, besides leading to a quantitative reduction 
of entanglement (compared to physical entanglement), 
demonstrates qualitative changes such as the possibil- 
ity of sudden death of effective entanglement in situ- 
ations when physical entanglement lives for arbitrarily 
long times. Our approach may be extended to estimate 
the time evolution of quantities describing quantum pro- 
cesses (e. g., entangling power or fidelity of a quantum 
process). 



II. THE MODEL 

The system under study consists of two parts, the mea- 
sured and the measuring subsystems. The former con- 
sists of the DQD ensemble, in which superpositions of 
excitonic states undergo pure dephasing due to the inter- 
action with the phonon modes of the surrounding crys- 
tal. The pure dephasing is only partial, as has been ex- 
perimentally shown [2H and later explained theoretically 
[Tgl . [20j : experimental and theoretical data yield both 
qualitative and quantitative agreement [20| . The |0) and 
|1) states of the qubit correspond to an empty QD and 
a QD with an exciton in its ground state, respectively 
(the two qubits are located in separate QDs). A num- 
ber of SETs which determine the charge distribution in 
their vicinity constitute the latter. This allows for the 
measurement of a set of elements of the DQD density 
matrix (the precise elements measured are determined 
by the geometry of the measurement device with respect 
to the DQD), and the subsequent calculation of effective 
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entanglement. 

The Hamiltonian which governs the evohition of the 
excitonic states in the DQD is 



H = ei(|l)(l|®I)+e2(I® |1)(1|) 

k 



(1) 



where the two states of each QD are denoted by |0) and 
|1), I is the unit operator, ei,2 are the transition energies 

(12) . 

in the two QDs, ' are exciton-phonon couphng con- 
stants, bl^jbk are creation and annihilation operators of 
the phonon modes, and uJk are the corresponding ener- 
gies (we put h = I). The explicit tensor notation refers 
to the DQD but is suppressed for the phonon reservoir 
components. 

Exciton wave functions are modelled by anisotropic 
Gaussians with the extension Zg/h in the xy plane for the 
electron/hole, and Iz along z for both particles. Then, 
the coupling constants for the deformation potential cou- 
pling between confined charges and longitudinal phonon 
modes have the form /i^'^'' = /fce*'*^^'^/^, where 



fk — 



2qVc 



V is the normalisation volume of the bosonic reservoir, 
d is the distance between the subsystems, k±^z are mo- 
mentum components in the xy plane and along the z axis, 
cre,h are deformation potential constants for electrons and 
holes, c is the speed of longitudinal sound, and g is the 
crystal density. 

In our calculations we use parameters typical for two 
self-assembled GaAs/InGaAs QDs stacked on top of each 
other [23, . The material parameters used are ctc = 
8 eV, (Th = —1 eV, c = 5.1 nm/ps, g = 5360 kg/m'^ 
(corresponding to GaAs), and d ~ 6 nm, 1^ = 4.4 nm, 
Zh = 3.6 nm, Iz — I nm (dot related parameters). 

The Hamiltonian can be diagonalised exactly using the 
Weyl operator method Q and the evolution is calculated 
following Ref. The interaction with the phonon 

modes leads to partial pure dephasing [13) HH, Ull , leaving 
the state occupations unchanged; the explicit forms of 
the time dependence of the off-diagonal density matrix 
elements may be found in Ref. [22]. 

The measuring subsystem is taken into account only in 
principle, in the sense that the information which can be 
gained about the state of the DQD ensemble is limited 
by realistic measurement capability. When measuring 
the state of a DQD by observing the current through a 
SET, the actual observable depends on its position with 
respect to the DQD. The number of the DQD density 
matrix elements which can be obtained is very limited; 
this restriction is a key point in this paper. There are a 
number of features inherent to this measurement scheme. 



e.g. finite measurement time and an uncertainty of the 
outcome, which are not taken into account for the sake 
of clarity. 

Let us consider a SET located near the lower QD (con- 
figuration A) . The interaction between the SET electron 
and the exciton in the dot shifts the energy levels in the 
SET when the QD is occupied, hence affecting the cur- 
rent. An appropriate choice of SET parameters allows for 
the maximisation of the difference in current flow and the 
measurement of the occupation of the lower dot [l^] • The 
measurement projectors corresponding to this situation 
are P„ = |0)(0| ® I and Pe = |1)(1| I, so the measured 
quantity is 



X = (00|p|00) + (01|p|01) 



(2) 



If such a SET is located near the upper QD, the measured 
quantity is 



y = (00|p|00)-t-(10|p|10). 



(3) 



Configuration B involves a SET located symmetrically 
between the QDs and in close proximity to them in such 
a way that the energy level on the SET island is sensitive 
to the probability of finding an exciton in midpoint TT] . 
The corresponding projectors are Pn = |— )(— | and Pe = 
|+)(+|, where 



|±) = (|01)±|10))/^/2. 



(4) 



This allows for the measurement of a linear combination 
of density matrix elements 



z = (01|p|01) + (10|p|10) + 2Re(01|p|10). 



(5) 



If the SET is located further away from the DQD re- 
gion (symmetrically), it is sensitive to the total number 
of excitons (configuration C). This allows for the mea- 
surement of 



d={n\p\u) 



(6) 



via projectors Pe = |11)(11| and P„ = |+)(-f | + |-)(-| + 
|00)(00|, or a = (00|p|00) via Pe = |00)(00| and P„ = 
l+>(+l + |->(-| + |ll)(ll| depending on the SET param- 
eters. Switching between the two modes can be accom- 
plished by applying different voltage to the SET. 



III. MINIMISING ENTANGLEMENT FROM 
THE SET DATA 

Consider the general quantum state p. There is a sim- 
ple 

Lemma .- For any convex entanglement measure E (or 
entanglement parameter) which is invariant under com- 
plex conjugate one has E{p) > i5(Re(p)), where the "Re" 
symbol means real part of the quantum state. 

Proof.- EiReip)) = E{^p+^p*) < ^E{p) + \E{p*) = 
E{p). 
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Any reasonable entanglement measure should be in- 
variant under the complex conjugate; all known entan- 
glement measures like entanglement of formation, con- 
currence, lognegativity, relative entropy of entanglement 
and all distillable quantities fulfil this condition. Some of 
them are also convex i.e. entanglement of formation, con- 
currence and relative entropy of entanglement. In what 
follows we shall consider the concurrence [26,, .27], which 
is both convex and invariant under the complex conju- 
gate. The reason for this choice of entanglement measure 
is its mathematical simplicity and the fact that it is eas- 
ily converted into Entanglement of Formation (which has 
a good physical interpretation). Since the SET measure- 
ment provides only real parts of off-diagonal density ma- 
trix elements, the first step of minimisation is to consider 
only the real parts of DQD states, p = Re(p). 

IV. SUDDEN DEATH OF EFFECTIVE 
ENTANGLEMENT 

Firstly, we will consider the situation, when a num- 
ber of SETs provide all possible information that can be 
gained about the state of the DQD with this measure- 
ment technique. This requires a pair of SETs in con- 
figuration A located near the two QDs, giving x and y 
defined, respectively, by Eq. ([2]) and dSj, and one in con- 
figuration C providing d (Eq. ©). Hence, all of the 
diagonal elements of the density matrix can be found. 



(OOIpIOO) 


= x + v 


+ d- 


\ — a 


(01|p|01> 


= i-y 


-d = 


b, 


(10|p|10> 


= l-x 


-d = 


c, 


(iiHii) 


^ d. 







A SET in configuration B provides z (Eq. ([5])) giving 
the real part of one of the off-diagonal density matrix 
elements 

Re((01|p|10)) = + 1 = Re(/i). (8) 

Since all of the diagonal density matrix elements are 
known, the set of initial maximally entangled states 
which cannot exhibit sudden death of effective entan- 
glement is the same as the set of states with real off- 
diagonal density matrix elements which do not exhibit 
sudden death of physical entanglement H^j- The time 
evolution of effective and physical entanglement in these 
states under pure dephasing is the same and the con- 
currence is equal to C{p) = 2Re(/i). The situation is 
different for states where all diagonal density matrix ele- 
ments are non-zero where the set of effectively entangled 
coherent states is substantially reduced. In Fig. [1] the 
effective entanglement of coherent states (Re(/i) — Vbc) 
as a function of b and c is plotted for two values of a. 

Secondly, let us consider the situation when only x 
and z have been measured (configuration A and B), so 
only linear combinations of some density matrix elements 




FIG. 1: The effective entanglement (colorscale) of coherent 
states in the first measurement setup as a function of known 
occupations b and c for two values of a. The plot shows only 
nonzero effective entanglement. 



are known. The measurement outcome x € [0,1]. The 
outcome z G [0,2], but it is easy to show that for non- 
zero effective entanglement z > I. The minimisation 
of entanglement requires that the quantity ad be max- 
imized, which occurs for Re(/i) ~ Vbc (the unknown 
density matrix elements a, 6, c, d — 1 — a — fe — c 
and Re{h) are defined in Eq. (O and ([5])); then ad = 
{x — b){l — X — {—Vb + -y/z)^). Finding the maximal ad, 
which leads to set values for a, b, c, d=l — a — b — c and 
Re(/i), and minimising over the five still unknown off- 
diagonal density matrix elements produces the effective 
entanglement for given measurement outcomes x and z. 
Effective entanglement as a function of x and z is plotted 
in Fig. H 

It is interesting to consider here the time-evolution of 
effective entanglement under phonon-induced pure de- 
phasing of an initially maximally entangled state |-|-) (Eq. 
(HI). When measured it will yield a; = 0.5 and z = 2, so 
the effective concurrence Ce(|-l-)(-f|) = 1 (note that the 
state |-) has x = 0.5, but z = and Ce(l-)(-|) = 0). 
Phonon-induced evolution of the state does not change x, 
but z decreases with decreasing Re(/i) leading to sudden 
death of effective entanglement for sufficiently dephased 
states. This state does not exhibit sudden death of phys- 
ical entanglement. 

Thirdly, we will consider the simplest situation, where 
the measurement device is limited to configuration B and 
the measurement data yield only z as defined in Eq. ([5]) . 
The amount of information gained by the measurement is 
very limited. Minimising effective entanglement requires 
the maximisation of ad the same as in the previous mea- 
surement setup, yet now the maximum is easily found 
for 6 = c = Re(/i) = z/4. The dependence of effective 
entanglement on z is plotted in the inset of Fig. [31 

The evolution of effective entanglement of the initial 
state |-|-) in this setup under realistic phonon-induced 
pure dephasing is plotted in Fig. [3] for different temper- 
atures. As is to be expected, effective disentanglement 
occurs faster than physical disentanglement. Further- 
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FIG. 2: The effective entanglement (colorscale) in the second 
measurement setup as a function of the observables x and z. 
The plot shows only nonzero effective entanglement. 




FIG. 3: The time-evolution of effective entanglement for dif- 
ferent temperatures. Inset: Effective entanglement as a func- 
tion of the measured parameter z 



V. CONCLUSION AND DISCUSSION 

We have considered the evolution of entanglement in 
a DQD system under phonon-induced pure dephasing 
minimized over a set of attainable data (the considered 
measurement schemes are based on SETs). In all of the 
considered measurement setups a reduction of effective 
(minimized) entanglement compared to the physical one 
was observed. Furthermore, in the setups where the 
knovifledge of QD occupations is limited, sudden death 
was found in situations when the physical entanglement 
lives for arbitrarily long times and its sudden death is 
not possible. For a limited range of temperatures (cou- 
pling strengths) even the sudden birth of entanglement 
occurs (due to a mechanism that does not cause sudden 
birth of physical entanglement). Hence, we have shown 
that the analysis of entanglement dynamics in systems 
with difficult measurement access can lead to qualitative 
differences in measured entanglement. 

The present analysis may provide a starting point to 
analogous research in general dynamics including quan- 
tum process tomography ([2^, Hll). For instance, one 
may want to probe the entangling power (30l | of a se- 
quence of quantum gates in time in a possibly cheap way, 
when only some of measurements are not costly. Then 
the principles of the above analysis may be applied with 
the help of the Choi- Jamiolkowski isomorphism [sil [s^j ■ 
Fundamental parameters of the dynamics, including the 
fidelity of a quantum process, can be then estimated, 
i.e. by minimisation in time under a restricted set of 
data, especially when coarse-graining information about 
the dynamics is needed. 



more, sudden death of entanglement appears for suffi- 
ciently high temperatures (e. g., when the dephasing is 
strong enough). For a limited range of temperatures, 
sudden birth of entanglement is also observed. The sec- 
ond phenomenon is due to the enhancement of coherence 
which occurs when wavepackets from the two QDs meet 
due to positive interference between them; this mecha- 
nism does not lead to the sudden birth of physical entan- 
glement [22j. 
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